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Abstract 


A direct search for lepton-flavour-violating (LFV) H —» /it decays of the recently dis¬ 
covered Higgs boson with the ATLAS detector at the LHC is presented. The analysis is 
performed in the H —» pr ha d channel, where Th a d is a hadronically decaying r-lepton. The 
search is based on the data sample of proton-proton collisions collected by the ATLAS ex¬ 
periment corresponding to an integrated luminosity of 20.3 fb -1 at a centre-of-mass energy 
of yfs = 8 TeV. No statistically significant excess of data over the predicted background 
is observed. The observed (expected) 95% confidence-level upper limit on the branching 
fraction, Br (H -» /it), is 1.85% (1.24%). 
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1 Introduction 


The observation of the Higgs boson [1,2] with a mass of about 125 GeV [3] by the ATLAS and CMS 
experiments is a great success of the Large Hadron Collider (LHC) physics program at CERN. The next 
important step in this program are searches for signs of new physics beyond the Standard Model (SM) 
and detailed studies of the Higgs boson properties. Direct evidence for physics beyond the SM could be 
indicated via lepton-flavour-violating (LFV) Higgs boson decays. If the SM is replaced with an effective 
field theory, which has a single Higgs boson and is required to be renormalizable only to a finite mass 
scale, then LFV couplings may be introduced [4]. LFV decays can also occur naturally in models with 
more than one Higgs doublet [5-8], composite Higgs models [9,10], models with flavour symmetries [11], 
Randall-Sundrum models [12] and many others [13-18]. 

There are three possibilities for LFV effects mediated via virtual Higgs bosons: n~e, r-yu, and x-e 
transitions. Indirect experimental constraints are reviewed and translated into constraints on Br(// —» 
e/u,iiT, ex) in recent papers [4,19]. Searches for /x —> ey [20] place a very stringent constraint on H —> e/x 
decays: Br (H —> e/x) < O(10“ 8 ) [4, 19]. The indirect constraints on H —> pr ,ex decays mostly come 
from searches for r —> /xy, ey [21-23] or other rare r-lepton decays [24], as well as from measure¬ 
ments of the anomalous magnetic moment of the muon and the electron [25] and are much less stringent: 
Br {H —> /ir, ex) < 0(10%) [4,19], A relatively large Br {H —> jxx) can be achieved without any particular 
tuning of the effective couplings, while a large Br (H —> ex) is possible only at the cost of some fine- 
tuning of the corresponding couplings [19]. It is also important to note that the presence of a H —> jxx 
signal would essentially exclude the presence of a H —> ex signal, and vice versa, at an experimentally 
observable level at the LHC due to strong experimental bounds on /x —> ey decays [19]. The CMS Collab¬ 
oration has recently performed the first direct search for LFV H —> yur decays [26] and reported a slight 
excess (2.4 standard deviations) of data over the predicted background. Their - results give a 1.57% upper 
limit on Br (H —> fxx) at the 95% confidence level (CL). 

This paper presents a search for lepton-flavour-violating H —* /xx decays of the recently discovered Higgs 
boson with one muon and one hadronically decaying r-lepton Chad) in the final state. The analysis is 
based on the data sample of pp collisions which was collected at a centre-of-mass energy of = § TeV 
and corresponds to an integrated luminosity of 20.3 I'fT 1 . 


2 The ATLAS detector and object reconstruction 

The ATLAS detector 1 is described in detail in ref. [27]. ATLAS consists of an inner tracking detector 
(ID) covering the range \tj\ < 2.5, surrounded by a superconducting solenoid providing a 2 T magnetic 
field, electromagnetic (|? 7 | < 3.2) and hadronic calorimeters (|;/| < 4.9) and a muon spectrometer (MS) 
(I 77 I < 2.7) with a toroidal magnetic field. 

The signature of LFV H —> /xx decays used in this search is characterised by the presence of an energetic 
muon, originating directly from a Higgs boson decay and carrying roughly half of its energy, and the had¬ 
ronic decay products of a r-lepton. The data were collected with a single-muon trigger with a transverse 


1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points 
upward. Cylindrical coordinates (r, (f >) are used in the transverse plane, 0 being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 6 as rj = — In tan(0/2). 
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momentum, py = p sin 6, threshold of p-y = 24 GeV. The H —> rr and the LFV H —>/./r signatures with 
a muon and Th a d in the final state share many common features. Therefore, the object definitions and 
data quality cuts used in this analysis are the same as those in the recently published ATLAS search for 
H —> tt decays [28]. A brief description of the object definitions is provided below. 

Muon candidates are reconstructed using an algorithm [29] that combines information from the ID and 
the MS. Muon quality criteria such as inner detector hit requirements are applied to achieve a precise 
measurement of the muon momentum and to reduce the misidentification rate. Muons are required to 
have p-y >10 GeV and to be within |//| < 2.5. Typical reconstruction and identification efficiencies for 
muons satisfying these selection criteria are above 95% [29]. Exactly one identified muon is required in 
this analysis. 

Electron candidates are reconstructed from energy clusters in the electromagnetic calorimeters matched 
to tracks in the ID. They are required to have a transverse energy, Ey = E sin 6, greater than 15 GeV, to 
be within the pseudorapidity range |//| < 2.47, and to satisfy the medium shower shape and track selection 
criteria defined in ref. [30]. Candidates found in the transition region between the end-cap and barrel 
calorimeters (1.37 < |//| < 1.52) are not considered. Typical reconstruction and identification efficiencies 
for electrons satisfying these selection criteria range between 80% and 90% depending on Ey and //. Since 
electrons do not appear in the H —» p + Th a d decay mode, events with identified electrons are rejected. 

Jets are reconstructed using the anfi-k t jet clustering algorithm [31] with a radius parameter R = 0.4, 
taking clusters of calorimeter cells with deposited energy as inputs. Fully calibrated jets [32] are required 
to be reconstructed in the range \q\ < 4.5 and to have p-y > 30 GeV. To reduce the contamination of jets 
by additional interactions in the same or neighbouring bunch crossings (pile-up), tracks originating from 
the primary vertex must contribute a large fraction of the jet p-y when summing the scalar p-y of all tracks 
in the jet. The primary vertex is chosen as the proton-proton collision vertex candidate with the highest 
sum of the squared transverse momenta of all associated tracks. This jet vertex fraction is required to be 
at least 50% for jets with \rj\ <2.4 and p-y <50 GeV (no cut is applied to jets with p-y > 50 GeV). Jets with 
no associated tracks are retained. In the pseudorapidity range \rj\ < 2.5, jets containing b-hadrons (A—jets) 
are selected using a tagging algorithm [33], which has an efficiency of ~70% for A—jets in tt events. The 
corresponding light-flavour jet misidentification probability is 0 . 1 - 1 %, depending on the p-y and 77 of the 
jet. Only a very small fraction of signal events have b- jets, therefore events with identified A-jets are 
vetoed in the selection of signal events. 

Hadronically decaying r-leptons are identified by means of a multivariate analysis technique [34] based 
on boosted decision trees, which exploits information about ID tracks and clusters in the electromagnetic 
and hadronic calorimeters. The Th a d candidates are required to have charge q T = ± 1 in units of electron 
charge, and must be 1- or 3-track (1- or 3-prong) candidates. Events with exactly one Th a d candidate 
satisfying the medium identification criteria [34] with py >20 GeV and |?/| <2.47 arc considered in 
this analysis. The identification efficiency for Th a d candidates satisfying these requirements is 55-60%. 
Dedicated criteria [34] to separate Th a d candidates from misidentified electrons are also applied, with a 
selection efficiency for true Th a d decays of 95%. To reduce the contamination due to backgrounds where 
a muon fakes a Th a d signature, events where an identified muon with p-y(p) > 4 GeV overlaps with an 
identified Th a d are rejected [28]. The probability to misidentify a jet with p-y > 20 GeV as a Th a d candidate 
is typically 1 - 2 %. 

The missing transverse momentum (with magnitude ii™ ss ) is reconstructed using the energy deposits in 
calorimeter cells calibrated according to the reconstructed physics objects (<?, y, Th a d, jets and p) with 
which they are associated [35]. The energy from calorimeter cells not associated with any physics object 
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is included in the Zs™ ss calculation. It is scaled by the ratio of the scalar sum of pj of tracks which 
originate from the primary vertex but are not matched to any objects and the scalar sum of p-\ of all 
tracks in the event which are not matched to objects. The scaling procedure achieves a more accurate 
reconstruction of £'™ 1SS in high pile-up conditions. In this search, Zf | mss is a signature of neutrinos, and it 
is used to select and reconstruct signal events, as described below. 


3 Event selection and categorization 

Signal H —» pr events in the pr h a d final state are characterised by the presence of an energetic muon and 
a Thad of opposite charge as well as moderate E™ ss , which tends to be aligned with the Th a d direction. 
Backgrounds for this signature can be broadly classified into two major categories: 

• Events with true muon and Thad signatures, dominated by irreducible Z/y* —* tt production with 
some contributions from the VV pr + X (where V = W,Z), tt, single-top and SM H —> tt 
production processes; these events exhibit a very strong charge correlation between the muon and 
the Thad; therefore, the expected number of events with opposite-sign (OS) charges (Nos) is much 
larger than the number of events with same-sign (SS) charges (Nss)- 

• Events with a fake Thad signature, dominated by VV'+jets events with some contribution from multi¬ 
jet (many multi-jet background events have genuine muons from semileptonic decays of heavy- 
flavour hadrons), diboson (VV), tt and single-top events with some charge asymmetry Nos > Nss', 
Z —> pp+ jets events, where a Thad signature can be faked by either a jet (no charge correlation) or 
a muon (strong charge correlation), also contribute to this category. 

Events with a fake Thad tend to have a much softer pT('nmd) spectrum and a larger angular separation 
between the Thad and £™ ss directions. These properties are exploited to suppress such backgrounds and 
define signal and control regions. Events with exactly one muon and exactly one Thad with pj(p) > 
26 GeV, px(Thad) > 45 GeV and \i](p) - ?7(Thad)l < 2 form a baseline sample. The \rj(p) - ?7(Thad)l cut 
has ~99% efficiency for signal and rejects a considerable fraction of multi-jet and VV'+jcts events. At this 
stage of the event selection, the identified muon is also required to be isolated [28] in the calorimeters 
and in the tracking detector in order to reduce contamination from the multi-jet background. Two signal 
regions are defined using the transverse mass, mj, 2 of the p-E™ lss and Th a d-£™ ss systems: OS events 
with mj(p, Ej lss ) > 40 GeV and n?x(Th a d, £'™ ss ) < 30 GeV form the signal region-1 (SRI), while OS 
events with mj(p, E'™ ss ) < 40 GeV and m-K+had, £'™ ss ) < 60 GeV form the signal region-2 (SR2). The 
dominant backgrounds in SRI and SR2 are VV'+jcts and Z/y* —> tt events, respectively. The model¬ 
ling of the VV'+jcts background is checked in a dedicated control region (WCR) formed by events with 
mj(p, £™ 1SS ) > 60 GeV and mxlThad, £™ ss ) > 40 GeV. The choice of the mj cuts to define SRI, SR2 and 
WCR is motivated by correlations between nij(p, £. , ’ 1ISS ) and mx(Th a d> E™ 1SS ) in // —> pr signal and major 
background (W+jets and Z/y* —> tt) events, as illustrated in figure 1. No events with identified ff-jets 
are allowed in SRI, SR2 and WCR. The modelling of the tt and single-top backgrounds is checked in a 
dedicated control region (TCR), formed by events that satisfy the baseline selection and have at least two 
jets, with at least one being ff-tagged. Table 1 provides a summary of the event selection cuts used to 
define the signal and control regions. 
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(1 - cos A0), where f = yu, Thad and A(p is the azimuthal separation between the directions of the lepton pi or 


Thad) and Ej ISS vectors. 
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Figure 1: Two-dimensional distributions of the transverse mass of the yU-E™ ss system, £™ ss ), and that of the 
r had-E™ ss system, Hi T (Thad, E™ ss ), in simulated Z/y* — » rr (top left plot), W+jets (top right plot), H —♦ \it signal 
(bottom left plot) and data (bottom right plot) events. All events are required to have a well-identified muon and 
Thad of opposite charge with pT(fhad) > 20 GeV and p i (j-i) > 26 GeV. 


The LFV signal is searched for by performing a fit to the mass distribution in data, m^. MC , reconstructed 
from the observed muon, Th a d and El ™ ss objects by means of the Missing Mass Calculator [36] (MMC). 
Conceptually, the MMC is a more sophisticated version of the collinear approximation [37]. The main 
improvement comes from requiring that the relative orientations of the neutrinos and other r-lepton decay 
products are consistent with the mass and kinematics of a r-lepton decay. This is achieved by maximising 
a probability defined in the kinematically allowed phase space region. The MMC used in the H —> 
tt analysis [28] is modified to take into account that there is only one neutrino from a hadronic r- 
lepton decay in LFV H —> qr events. For a Higgs boson with mu = 125 GeV, the reconstructed m^. MC 
distribution has a roughly Gaussian shape with a full width at half maximum of ~19 GeV. The analysis is 
performed “blinded” in the 110 GeV< <150 GeV regions of SRI and SR2, which contains -94% 

of the expected signal events. 


4 Background estimation 


The background estimation method takes into account the background properties and composition dis¬ 
cussed in section 3. It also relies on the assumption that the shape of the m^ MC distribution for the 
multi-jet background is the same for OS and SS events. This assumption was verified in the published 
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Table 1: Summary of the event selection criteria used to define the signal and various control regions (see text). 


Cut 

SRI 

SR2 

WCR 

TCR 

Pt(m) 

>26 GeV 

>26 GeV 

>26 GeV 

>26 GeV 

PTOhad) 

>45 GeV 

>45 GeV 

>45 GeV 

>45 GeV 

m T (p,E™ ss ) 

>40 GeV 

<40 GeV 

>60 GeV 

- 

Wt(t had, -E™ 88 ) 

<30 GeV 

<60 GeV 

>40 GeV 

- 

Wd) ~ 7<Thad)l 

<2 

<2 

<2 

<2 

A'jet 

- 

- 

- 

>1 

N fo-jet 

0 

0 

0 

>0 


H —> tt search [38]. In addition, it was confirmed using a dedicated control region, MJCR, with an 
enhanced contribution from the multi-jet background. Events in this control region are required to pass 
all criteria for SRI and SR2 with the exception of the requirement on - //(rh a d)l, which is reversed: 
|- ;/(T| ia d)| > 2. Therefore, the number of the total OS background events, in each bin of the 
/?7 m mc ( or an y other) distribution in SRI and SR2 can be obtained according to the following formula: 


N 


bkg 


OS 


= fQCD 




TT 

-ss 


+ v 




os- 


N W+jets „top 
SS ^^OS-SS ^ JV OS- 


ss 


+ N, 


vv 

os-ss 


+ Ni 


•TT 

-SS’ 


(1) 


where the individual terms are described below. is the number of SS data events, which are domin¬ 
ated by VV'+jets events but also contain contributions from multi-jet and other backgrounds. The fractions 
of multi-jet background in SS data events inside the 110 GeV< m^ MC <150 GeV mass window are 

~17% and ~44% in SRI and SR2, respectively. The contributions Aq^Lss = ^os § ' ~ r QCD ■ (Vgg g ' 
are add-on terms for the different backgrounds components (where bkg-/ indicates the / th background 
source: Z —> tt, Z —> /j/u, IV+jets, VV, H —> tt and events with f-quarks), which also account for com¬ 
ponents of these backgrounds already included in SS data events. 3 * * The factor tqqd - /V™' 11 ' 3et /Af™ ulti Jel 
accounts for potential differences in flavour composition (and, as a consequence, in jet —» Th a d fake 
rates) of final-state jets introduced by the same-sign or opposite-sign charge requirements. The value of 
rQCD = 1.10 ± 0.14 is obtained from a multi-jet-enriched control region in data, as discussed in detail in 
ref. [38]. It was verified in the MJCR in this search. 


The largely irreducible Z/y* —» tt background is modelled by Z/y* —> /iji data events, where the muon 
tracks and associated energy deposits in the calorimeters are replaced by the corresponding simulated 
signatures of the final-state particles of the r-lepton decay. In this approach, essential features such as 
the modelling of the kinematics of the produced boson, the modelling of the hadronic activity of the 
event (jets and underlying event) as well as contributions from pile-up are taken from data. Therefore, 
the dependence on the simulation is minimized and only the r-lepton decays and the detector response 
to the r-lepton decay products are based on simulation. This hybrid sample is referred to as embedded 
data in the following. A detailed description of the embedding procedure can be found in ref. [39]. The 
Z/y* —> tt normalization is a free-floating parameter in the final fit to data and it is mainly constrained 
by events with 60 GeV<m^. MC <l 10 GeV in SR2. 

The VT+jets and Z —> pp backgrounds arc modelled by the ALPGEN [40] event generator interfaced with 
PYTHIA8 [41] to provide the parton showering, hadronization and the modelling of the underlying event. 

3 The j-qcd • A^ss 8 ' correction in the add-on term is needed because same-sign data events include multi-jet as well as elec- 

troweak events (Z —> tt, Z — > /a/a, Vk+jets, VV, H —> tt and events with t-quarks) and their contributions cannot be 

separated. 
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In all VV'+jets events, the Th a d signature is faked by jets. The WCR region is used to check the modelling 
of the VV'+jcts kinematics and to obtain normalizations for OS and SS VV'+jcts events. An additional 
overall normalization factor for the /V ( ^^ term in eq. (1) is introduced as a free-floating parameter in 
the final fit in SRI. By studying WCR events and SRI events with m“ MC > 150 GeV (dominated by 
W +jcts background), it is also found that a m^ MC shape correction, which depends on the number of jets, 
/AOhad) and | rj(ju) - //(r[ ia d)l, needs to be applied in SRI. This correction is derived from SRI events with 
m M MC > 15 Q GeV and it is applied to events with all values of m“ MC . A 50% difference between the 
SRl-based correction and that obtained in WCR is taken as a corresponding modelling uncertainty on the 
m M MC s [ ia p C f or q lc IT +jcts background in SRI. The size of this uncertainty depends on mJ^ MC and it is 
as large as ±10% for VV'+jcts events with m“ MC < 150 GeV. In the case of SR2, a good modelling of the 
/Vjet, pxGhad) and | rj(ju) - ?7(Th a d)| distributions suggests that such a correction is not needed. However, a 
modelling uncertainty on the m^ MC shape of the VV'+jcts background in SR2 is assigned based on the 50% 
difference between the default n?^ MC shape and the one obtained after applying the correction derived for 
SRI events. The size of this uncertainty is below 5% in the 110 GeV< m^. MC <150 GeV region, that 
contains most of the signal events. It was also checked that applying the same correction in SR2 as in SRI 
had a negligible effect on the final result and the extracted branching ratio Br(// —> pr) (see Section 6) 
would only be affected at a level below 3%. The modelling of jet fragmentation and the underlying event 
has a significant effect on the estimate of the jet —» Th a d fake rate in different regions of the phase space and 
has to be accounted for with a corresponding systematic uncertainty. To estimate this effect, the analysis 
was repeated using a sample of VV'+jets events modelled by ALPGEN interfaced with the HERWIG [42] 
event generator. Differences in the VV'+jets predictions in SRI and SR2 are found to be ±9% and ±2%, 
respectively, and are taken as corresponding systematic uncertainties. 

In the case of the Z —» qp background, there are two components: events where a muon fakes a Th a d and 
events where a jet fakes a Th a d- Predictions for the shape and normalization of the Z —> /j/u (where /i —» 
T hlicT) background are obtained from simulation. For the Z —» fi/i background where a jet is misidentified 
as a Th a d candidate, the normalization factor and shape corrections, which depend on the number of jets, 
p-rGhad) and |//(p) - ;/(r| ia d)|, are derived by using events with two identified OS muons with an invariant 
mass, m Mfl , in the range of 80-100 GeV. Since this background does not have an OS-SS charge asymmetry, 
a single correction factor is derived for OS and SS events. A 50% difference between the mJ^ MC shape 
with and without this correction is taken as a corresponding systematic uncertainty. 

The backgrounds with top quarks arc modelled by the POWHEG [43-45] (for tt, Wt and ^-channel 
single-top production) and AcerMC [46] (/-channel single-top production) event generators interfaced 
with PYTHIA8 to provide the parton showering, hadronization and the modelling of the underlying event. 
The TCR is used to check the modelling and to obtain normalizations for OS and SS events with top 
quarks. The normalization factors obtained in the TCR are extrapolated into SRI and SR2, where tt 
and single-top events may have different properties. To estimate the uncertainty associated with such an 
extrapolation, the analysis is repeated using the MC@NLO [47] event generator instead of POWHEG for 
tt production. 4 This uncertainty is found to be ±7.2% (±3.7%) for backgrounds with top quarks in SRI 
(SR2). 

The background due to diboson {WW, ZZ and WZ) production is estimated from simulation, normalized 
to the cross sections calculated at next-to-leading order (NLO) in QCD [48]. The ALPGEN event gen¬ 
erator interfaced with HERWIG is used to model the WW process, and HERWIG is used for the ZZ and 
WZ processes. 


4 The same extrapolation uncertainty is assumed for tt and single-top backgrounds. 
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Figure 2: Distributions of the mass reconstructed by the Missing Mass Calculator, m^. MC , in SRI (left) and SR2 
(right). The background distributions are determined in a global fit. The signal distribution corresponds to Br (H —> 
fjT)-25%. The bottom panel of each sub-figure shows the ratio of the observed data and the estimated background. 
The grey band for the ratio illustrates post-fit systematic uncertainties on the background prediction. The statistical 
uncertainties for data and background predictions are added in quadrature for the ratios. The last bin in each 
distribution contains overflow events. 


Finally, events with Higgs bosons produced via gluon fusion or vector-boson fusion (VBF) processes are 
generated at NLO accuracy with POWHEG [49] event generator interfaced with PYTHIA8 to provide the 
parton showering, hadronization and the modelling of the underlying event. The associated production 
(ZH and WH) samples are simulated using PYTHIA8. All events with Higgs bosons are produced with a 
mass of mu = 125 GeV and normalized to cross sections calculated at next-to-next-to-leading order in 
QCD [50-52]. The SM H —» rr decays are simulated by PYTHIA8. The LFV Higgs boson decays are 
modelled by the EvtGen [53] event generator according to the phase-space model. In the H —» /ut decays, 
the r-lepton decays are treated as unpolarised because the left- and right-handed r-lepton polarisation 
states are produced at equal rates. 

All simulated samples are passed through the GEANT4-based ATLAS detector simulation [54,55]. The 
simulated events are overlaid with additional minimum-bias events to account for the effect of multiple 
pp interactions (pile-up) occurring in the same and neighbouring bunch crossings. 

Figure 2 shows the m^ MC distributions for data and the predicted backgrounds in each of the signal re¬ 
gions. The backgrounds are estimated using the method described above. The signal efficiencies for 
passing the SRI or SR2 selection requirements are 2.1% and 1.5%, respectively, and the combined ef¬ 
ficiency is 3.6%. The numbers of observed events in the data as well as the signal and background 
predictions in the mass region 110 GeV< m“ MC <150 GeV can be found in table 2. 
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Table 2: Data yields, signal and post-fit OS-SS background predictions (see eq. (1)) for the 110 GeV< 
m“ MC <150 GeV region. The signal predictions are given for Br (H —> /rr)=0.77%. The background predictions 
are obtained from the combined fit to SRI, SR2, WCR and TCR. The post-fit values of systematic uncertainties 
are provided for the background predictions. For the total background, all correlations between various sources of 
systematic uncertainties and backgrounds are taken into account. The quoted uncertainties represent the statistical 
(first) and systematic (second) uncertainties, respectively. 



SRI 

SR2 

Signal 

69.1 ±0.8 ±9.2 

48.5 ±0.8 ±7.5 

Z —> TT 

133.4 ±6.9 ±9.1 

262.6 ± 9.7 ± 18.6 

W +j ets 

619 ±54 ±55 

406 ± 42 ± 34 

Top 

39.5 ±5.3 ±4.7 

19.6 ±3.1 ±3.3 

Same-Sign events 

335 ± 19 ±47 

238 ± 16 ± 34 

VV + Z^nn 

90 ±21 ± 16 

81 ±22± 17 

H —» TT 

6.82 ±0.21 ±0.97 

13.7 ± 0.3 ± 1.9 

Total background 

1224 ± 62 ± 63 

1021 ±51 ±49 

Data 

1217 

1075 


5 Systematic uncertainties 

The largest systematic uncertainties arise from the normalization (±10% uncertainty) and modelling 5 of 
the VV'+jcts background. The uncertainties on rQCD (±12.7%) and on the normalization (±6% uncer¬ 
tainty) and modelling of Z —* rr also play an important role. The other major sources of experimental 
uncertainty, affecting both the shape and normalization of signal and backgrounds, are the uncertainty 
on the Thad energy scale [34] (measured with ±(2-4) % precision) and uncertainties on the embedding 
method used to model the Z —» rr background [28]. Less significant sources of experimental uncertainty, 
affecting the shape and normalization of signal and backgrounds, are the uncertainty on the jet energy 
scale [32,56] and resolution [57], The uncertainties in the Th a d energy resolution, the momentum scale 
and resolution of muons, and the scale uncertainty on E™ ss due to the energy in calorimeter cells not 
associated with physics objects are taken into account, however, they are found to be relatively small. 
The following experimental uncertainties primarily affect the normalization of signal and backgrounds: 
the ±2.8% uncertainty on the integrated luminosity [58], the uncertainty on the Thad identification effi¬ 
ciency [34], which is measured to be ±(2—3)% for 1-prong and ±(3—5)% for 3-prong decays, the ±2.1% 
uncertainty for triggering, reconstructing and identifying muons [29,59], and the ±2% uncertainty on the 
Z?—jet tagging efficiency [33], 

Theoretical uncertainties are estimated for the signal and for the H — > rr, VV and Z —> juju (with ju —» 
rS e ) backgrounds, which are modelled with the simulation and are not normalized to data in dedicated 
control regions. Uncertainties due to missing higher-order QCD corrections on the production cross 
sections are found to be [60] ±10.1% (±7.8%) for the Higgs boson production via gluon fusion in SRI 
(SR2), ±1% for the Z —» ju/r background and for VBF and VH Higgs boson production, and ±5% for 
the VV background. The systematic uncertainties due to the choice of parton distribution functions used 
in the simulation are evaluated based on the prescription described in ref. [60] and the following values 
are used in this analysis: ±7.5% for the Higgs boson production via gluon fusion, ±2.8% for the VBF 

5 Some of these uncertainties (e.g., uncertainties due to shape corrections and extrapolation uncertainties) are discussed 

in the text above. 
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Figure 3: Post-fit combined m^ MC distribution obtained by adding individual distributions in SRI and SR2. In 
the lower part of the figure, the data are shown after subtraction of the estimated backgrounds. The grey band 
in the bottom panel illustrates the post-fit systematic uncertainties on the background prediction. The statistical 
uncertainties for data and background predictions are added in quadrature on the bottom part of the figure. The 
signal is shown assuming Br(7/ —> pr)=0.77%, the central value of the best fit to Br (H —> /jt). The last bin of the 
distribution contains overflow events. 


and VH Higgs boson production, and ±4% for the Z —* p/r and VV backgrounds. Finally, an additional 
±5.7% systematic uncertainty on Br (H —> tt) is applied to the SM H —> rr background. 


6 Results 

A simultaneous binned maximum-likelihood fit is performed on the m“ MC distributions in SRI and SR2 
and on event yields in WCR and TCR to extract the LFV branching ratio Br (H —> iut). The fit exploits 
the control regions and the distinct shapes of the VT+jets and Z —» rr backgrounds in the signal regions to 
constrain some of the systematic uncertainties. This leads to an improved sensitivity of the analysis. The 
post-fit /n^ MC distributions in SRI and SR2 are shown in figure 2, and the combined /n“ MC distribution 
for both signal regions is presented in figure 3. Figure 2 illustrates good agreement between data and 
background expectations in SRI. A small excess of the data over the predicted background is observed 
in the 120 GeV< m^ MC <140 GeV region in SR2. This small excess in SR2 has a local significance of 
2.2 standard deviations and a combined significance for both signal regions of 1.3 standard deviations. 
This corresponds to a best fit value for the branching fraction of Br (H —> iut)=(0.11 ± 0.62)%. Due to 
the low significance of the observed excess, an upper limit on the LFV branching ratio \ii (H —> /it) for 
a Higgs boson with mu = 125 GeV is set using the CL V modified frequentist formalism [61] with the 
profile likelihood-ratio test statistics [62]. The observed and the median expected 95% CL upper limits 
are 1.85% and 1.24+jj^ 35 %, respectively. Table 3 provides a summary of all results. 
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Table 3: The expected and observed 95% confidence level (CL) upper limits and the best fit values for the branching 
fractions for the two signal regions and their combination. 



SRI 

SR2 

Combined 

Expected limit on Br (H —» /it) [%] 

1 60 +0 - 64 
-0.45 

1 75 +0 - 71 
• - 0.49 

I 24 +0 ' 50 

II -0.35 

Observed limit on Br (H —» qr) [%] 

1.55 

3.51 

1.85 

Best fit Br (H -» qr) [%] 

-0 07 +0 - 81 
ua - 0.86 

1 94+0.92 
L,y -0.S9 

0.77±0.62 


7 Summary 


A direct search for lepton-flavour-violating H —> fir decays of the recently discovered Higgs boson is 
performed in the Thad decay mode of r-leptons using a data sample of proton-proton collisions recorded 
by the ATLAS detector at the LHC corresponding to an integrated luminosity of 20.3 fb -1 at a centre- 
of-mass energy of yfs = 8 TcV. The observed and the median expected upper limits at 95% CL on the 
branching fraction, Br (H —» /it), are 1.85% and 1.24 1 [[ 35 %, respectively. This search places significantly 
more stringent constraints on Br {H —> pr) compared to earlier indirect estimates. The result of this 
analysis is found to be consistent with the one published by the CMS Collaboration [26]. 
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